Nucleon strangeness form factors and moments of PDF 
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Abstract. 

The calculation of the nucleon strangeness form factors from Nt ■ = 2 + 1 clover fermion lattice QCD is presented. 
Disconnected insertions are evaluated using the Z(4) stochastic method, along with unbiased subtractions from the hopping 
parameter expansion. We find that increasing the number of nucleon sources for each configuration improves the signal 
significantly. We obtain G S M (0) = —0.017(25) (07), which is consistent with experimental values, and has an order of 
magnitude smaller error. Preliminary results for the strangeness contribution to the second moment of the parton distribution 
function are also presented. 

PACS: 13.40.-f, 12.38.Gc, 14.20.Dh 



INTRODUCTION 

Understanding the structure of the nucleon from QCD has been one of the central issues in hadron physics. In 
particular, the strangeness content of the nucleon attracts a great deal of interest lately. It is also an ideal probe for 
the virtual sea quarks in the nucleon. Extensive experimental/theoretical studies indicate that the strangeness content 
varies depending on the quantum number carried by the ss pair: the scalar density is about 0-20% of that of up, 
down quarks, the quark spin is about —10 to 0% of the nucleon, and the momentum fraction is only a few percent of 
the nucleon. In general, the uncertainties in the strangeness matrix elements are quite large in both experiments and 
theories. Under these circumstances, it is desirable to provide the definitive quantitative results using lattice QCD. 

The challenge in the lattice QCD calculation of strangeness matrix elements resides in the evaluation of the so-called 
disconnected insertion (DI). In fact, it requires the calculation of all-to-all propagators, which is prohibitively expensive 
compared to the connected insertion (CI). Consequently, there are only a few DI calculations [1, 2, 3], where the all- 
to-all propagators are stochastically estimated [4]. In this proceeding, we report the improvement of the calculation 
of all-to-all propagators using the stochastic method along with unbiased subtractions from the hopping parameter 
expansion [5], and the increment of the number of nucleon sources [6, 7]. We present the results for the strangeness 
contribution to the electromagnetic form factors [7] and the second moment of the nucleon. The preliminary result for 
the first moment of the nucleon is presented in Ref. [8]. 

FORMALISM AND SIMULATION PARAMETERS 

We employ Nf = 2+1 dynamical configurations with nonperturbatively &{a) improved clover fermion and RG- 
improved gauge action generated by CP-PACS/JLQCD Collaborations [9]. We use /3 = 1.83 and c sw = 1.7610 
configurations with the lattice size of I? x T = 16 3 x 32, which corresponds to (2fm) 3 box in physical spacial size 
with the lattice spacing of cr l = 1.625GeV [9]. For the hopping parameters of u,d quarks (JQ^) and s quark (k s ), we 
use K ud = 0.13825, 0.13800, and 0.13760, which correspond to m n = 0.60, 0.70, and 0.84 GeV, respectively, and 
K s = 0.13760 is fixed. We perform the calculation only at the dynamical quark mass points, where 800 configurations 
are used for K ud = 0.13760, and 810 configurations for K ud = 0.13800, 0.13825. 



The nucleon matrix elements can be obtained through the calculation of 3pt function Tlf (as well as 2pt function 
n 2 P'), defined by 

n 3 /\p,t 2 ; q, h ; p' = p-q,t ) = £ . e +*&-*>) (0\T [XN&,t 2 )J(ZiA)XN$o,to)} |0), (1) 

where %n is the nucleon interpolating field and / is the insertion operator. Since there is no strange quark as a valence 
quark in the nucleon, the 3pt is a DI which entails a multiplication of the nucleon 2pt correlator with the current quark 
loop. For the evaluation of the quark loop, we use the stochastic method [4], with Z(4) noises in color, spin and space- 
time indices. We generate independent noises for different configurations, in order to avoid possible auto-correlation. 
We use Nnoise — 600 noises for k uc \ — 0.13760,0.13800 and N m ,i se — 800 for k u( { — 0.13825. To reduce fluctuations, 
the charge conjugation and 75-hermiticity (CH), and parity symmetry are used [6, 7]. We also perform unbiased 
subtractions [5] to reduce the off-diagonal contaminations to the variance. For subtraction operators, we employ those 

obtained through hopping parameter expansion (HPE) for the propagator M~ 1 , j^M~ l = + j^{kD)j^ H 

where D denotes the Wilson-Dirac operator and C the clover term. We subtract up to order (kD) a ((fcD) 3 ) term for the 
form factor (second moment) calculation, and observe that the statistical errors become about 50 (70) %, compare to 
the results without subtraction. 

In the stochastic method, it is quite expensive to achieve a good signal to noise ratio (S/N) just by increasing N n0 i Se 
because S/N improves with \/N no j se . In view of this, we use many nucleon point sources N src in the evaluation of the 
2pt part for each configuration. Since the calculations of the loop part and 2pt part are independent of each other, this is 
expected to be an efficient way. We take N src = 64 for K ud = Q.\ 3760 and N src = 82 for K Ild — 0. 1 3800, 0. 1 3825, where 
locations of sources are taken so that they are separated in 4D-volume as much as possible. Details of the simulation 
setup are given in Ref . [7] . 



STRANGENESS ELECTROMAGNETIC FORM FACTORS 

The formulas for Sachs electric (magnetic) form factors G S E {G S M ) are given by 
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where J fl (x + ji/2) = j q{x)(l — y^)U il {x)q{x + ji) — q(x + \ + Yii)U^(x)q(x) is the point-split conserved vector 



operator, {ij,k} ^ 4, = (1 ± 74)/2 , = (±i)/2 x (1 ±74)757^ and£^ = ^Jm^ + q 2 . The upper sign coiTesponds 
to the forward propagation (f 2 ^> h ^ to), and the lower sign corresponds to the backward propagation (f 2 <C f 1 to). 
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FIGURE 1. R' M (left) and R' E (right) with K ud = 0. 13760, N src = 64 (circles) and N src = 4 (triangles), plotted against the nucleon 
sink time f 2 - The dashed line is the linear fit where the slope corresponds to the form factor. 
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FIGURE 2. The chiral extrapolated results for G S M (Q 2 ) (left) and G S E (Q 2 ) (right) plotted with solid lines. Shaded regions 
represent the statistical and systematic error added in quadrature. Shown together are the lattice data for each K uc j. 



In Fig. 1 , we plot typical figures for R' M E , where R' M E = —J— Y!u=L+t w i m E being trivial kinematic factors 

' K M,E * 

in Eq. (3). Since R 1 ME = const. + ?2 x G i ME , the linear slope corresponds to the signal of the form factor. One can 
observe the significant S/N improvement by increasing N src - In fact, the improvement is found to be nearly a factor of 
\fN S rc (ideal improvement). 

We then study the Q 2 dependence of the form factors. For the magnetic form factor, we employ the dipole form, 
G'm(Q 2 ) — G S M (Q)/(1 +Q 2 /A 2 ) 2 , where reasonable agreement with lattice data is observed. For the electric form 
factor, we employ G E (Q 2 ) = g s E ■ Q 2 /(l +Q 2 /A 2 ) 2 , considering that G E (0) = from the vector current conservation, 
and the pole mass A is taken from the fit of magnetic form factor. 

Finally, we perform the chiral extrapolation for the fitted parameters. Since our quark masses are relatively heavy, 
we consider only the leading dependence on trig, which is obtained by heavy baryon chiral perturbation theory 

(HB^PT) [10]. The chiral extrapolated results are G S M (0) = -0.017(25), Aa = 0.58(16), (rj) M = -6^p-| 2 2 =0 = 
-7.4(71) x 10- 3 fm 2 and g| =0.027(16) (or (r 2 ) E = -6^\ Q 2 =0 = -2.4(15) x 10- 3 fm 2 ). 

We examine the systematic uncertainties in the result of form factors. For the ambiguity of Q 2 dependence, we 
reanalyze the data using the monopole form, and obtain the results which are consistent with those from the dipole 
form. For the uncertainties in chiral extrapolation, we test two alternative extrapolations [7], and find that all results 
are consistent with each other. For the contamination from excited states, we employ the new projection operator [7] 
which eliminates the Su state, and conclude that such contaminations are negligible. 

Our final result for the magnetic moment is G^(0) = —0.017(25) (07), where the first error is statistical and 
the second is systematic from uncertainties of the Q 2 extrapolation and chiral extrapolation. We also obtain Aa = 
0.58(16)(19) for dipole mass or Aa = 0.34(17)(11) for monopole mass, and g% = 0.027(16)(08). These lead to 
Gm(Q 2 ) = -0.015(23), G E (Q 2 ) =0.0022(19) at Q 2 =0.1GeV 2 , where error is obtained by quadrature from statistical 
and systematic errors. In Fig. 2, we plot G S M (Q 2 ), G E (Q 2 ), where the shaded regions correspond to the square-summed 
error. Compared to the global analysis of the experimental data, e.g., G S M (Q 2 ) = 0.29(21) and G S E (Q 2 ) = —0.008(16) 
at2 2 = 0.1GeV 2 [11], our results are consistent with them, with an order of magnitude smaller error [7] . 



SECOND MOMENT OF THE NUCLEON 

The (asymmetry of) strangeness second moment of the nucleon (x 2 ) iS _f = J Q dxx 2 {s(x) — s(x)) can be obtained by 
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FIGURE 3. LEFT: The ratio of 3pt to 2pt with K ud = 0.13760, N src = 64 (circles) and N src = 4 (triangles), plotted against the 
nucleon sink time ti. The dashed line is the linear fit where the slope corresponds to the second moment. RIGHT: The lattice bare 
results for the second moment at each valence quark mass K uc / for the nucleon, plotted against (mxa) 2 . The dashed line corresponds 
to the linear chiral extrapolation, and the red point is the chiral extrapolated result. 



with the three-index operator defined as 
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where i ^ 4, and the upper (lower) sign corresponds to the forward (backward) propagation as before. 

In Fig. 3 (left), we plot the ratio of 3pt to 2pt for (jc 2 ). s _.v in terms of t% for k uc i = 0.13760, p 2 = {2n/La) 2 , where 
the summation of operator insertion time t\ is taken as was done for the form factor analysis. Note that the linear slope 
corresponds to the signal for (x 2 ) s -s- One can clearly see that increasing N src reduces the error bar significantly again 
(about a factor of y/N src , i.e., almost ideally). In Fig. 3 (right), we plot the bare value of the (x 2 ) 4 _ s - in terms of (mxa) 2 , 
and perform the chiral extrapolation. We find that the result at each K M / and the chiral extrapolated result are basically 
consistent with zero within the error-bar. For the final quantitative result, it is necessary to take the renormalization 
factor into account. Systematic uncertainties have to be examined as well. The study along this line is in progress. 
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